Testing the process dependence of the Sivers function 
via hadron distributions inside a jet 
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We study the process dependence of the Sivers function by considering the impact of color-gauge 
invariant initial and final state interactions on transverse spin asymmetries in proton-proton scat- 
tering reactions within the framework of the transverse momentum dependent (TMD), generalized 
parton model. To this aim, we consider the azimuthal distribution of leading pions inside a fragment- 
ing jet as well as single inclusive jet asymmetry in polarized proton-proton collisions. In contrast 
to single inclusive pion production, in both cases we can isolate the Sivers contribution and thereby 
study its process dependence. The predictions for the Sivers asymmetry obtained with and without 
inclusion of color gauge factors are comparable in size but with opposite signs. We conclude that 
both processes represent unique opportunities to discriminate among the two approaches and test 
the universality properties of the Sivers function in hadronic scattering reactions. 

PACS numbers: 12.38.Bx, 13.85.Ni, 13.88.+C 



Single transverse-spin asymmetries (SSAs) in high en- 
ergy lepton-hadron and hadronic scattering processes 
have garnered considerable attention from both exper- 
imental and theoretical communities [l[ . Generally, they 
are denned as the ratio of the difference and the sum of 
the cross sections when the hadron's spin vector S± is 
flipped, A N = (a(S ± ) - a(-S 1 _))/(a(S±) + <r(-S±)) = 
Acr/(2<7 unp ). The SSAs for single inclusive particle pro- 
duction in proton-proton scattering are among the earli- 
est processes studied Q and remain extremely challeng- 
ing to explain in the context of perturbative quantum 
chromodynamics (QCD) 3]. The trend of large SSAs 
in the pioneering fixed target experiments has been ob- 
served over a wide range of energies and more recently at 
significantly larger center-of-mass energies in the proton- 
proton collision experiments at Relativistic Heavy Ion 
Collider (RHIC) [4j, |5j. Also, azimuthal and transverse- 
spin asymmetries have been observed in Drell-Yan (DY) 
processes IfiL m semi-inclusive deep inelastic scattering 
(SIDIS) [3j8(j and in hadron pair production in e + e~ 
scattering [9|. 

From a theoretical perspective SSAs are characterized 
by the interference between helicity flip and non-flip scat- 
tering amplitudes with a relative color phase. Two ap- 
proaches have been proposed in the framework of per- 
turbative QCD to account for these effects. On the one 
hand is the collinear factorization formalism at next-to- 
lcading-power (twist-3) in the hard scale where SSAs 
are given by a convolution of universal non-perturbative 
quark-gluon-quark correlation functions and hard scat- 
tering amplitudes [loWl^. 



tation functions (PDFs and FFs). Prominent examples 
are the quark Sivers function [1 31 ] . which represents the 
azimuthal distribution of unpolarized quarks in a trans- 
versely polarized nucleon and the Collins fragmentation 
function [Tij ]. which describes the production of pseudo- 
scalar mesons (or unpolarized hadrons) from transversely 
polarized fragmenting quarks. In this approach color 
phases are given by initial and/or final state interactions 
(ISIs/FSIs) between the active quark and spectator rem- 
nants in the full scattering amplitude. The details of the 
ISIs and FSIs depend on the scattering process and for 
PDFs such as the Sivers function, these color phases are 
incorporated into the Wilson lines of the gauge invariant 
definition of TMD PDFs. It is a fundamental prediction 
of QCD factorization that the form of the gauge link de- 
pends on the hard sub-process 15] indicating that the 
Sivers function is non-universal [16J. The oft-discussed 
case is the difference between the FSIs in SIDIS and the 
ISIs in DY scattering which leads to the prediction of 
an opposite relative color factor [la ]. Further, applying 
similar reasoning to hadron production in proton-proton 
collisions, typically the Sivers function has a more com- 
plicated color factor structure since both ISIs and FSIs 
contribute flBT-tlil, 



The other framework relies on factorization in terms 
of a hard scattering cross section and transverse momen- 
tum dependent (TMD) parton distribution and fragmen- 



While TMD factorization has not been established for 
hadron production in hadronic reactions 19], an exten- 
sive program of phenomenology has been carried out by 
including the correlations of intrinsic parton motion and 
transverse spin in the context of the so-called generalized 
parton model (GPM) . Introduced [2(j as a generalization 
of the collinear perturbative QCD approach, it has been 
used to describe the SSAs for inclusive particle produc- 
tion 2l|. Here factorization has been assumed as a rea- 
sonable starting point for analyses. At the same time, 
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the leading-twist naive time-reversal odd (T-odd) TMD 
PDFs have conditionally been assumed to be universal. 

In this Letter we present an analysis of SSAs in proton- 
proton scattering while taking into account the effects of 
ISIs and FSIs and allowing for process dependence within 
the framework of GPM. This will be referred to as the 
color gauge invariant (CGI) GPM. Previous studies along 
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these lines have been carried out in 

We concentrate on reactions where one can explore the 
crucial issue of process dependence and universality of the 
Sivers function. Since several competing mechanisms can 
play a role in hadron collisions, following [22| we consider 
the process p^p — > jet tt + X , where one observes a large 
Pt jet and looks for the azimuthal distribution of leading 
pions within the jet: a process under active investigation 
by the STAR Collaboration at RHIC [H]. By contrast, 

I 



an analysis where transverse partonic motion was consid- 
ered only in the fragmentation process, aimed at a study 
of the universality of the Collins effect, was presented 



24j. It is also important to note that the process we 



are studying is different from the case where two almost 
back- to-back hadrons or jets are observed as in Til 17|. 
A single jet is measured in our study. Thus, the a naly sis 
of process dependence follows that carried out in |18j |. 



Compared with inclusive pion production, we empha- 
size that the Sivers and Collins contributions can be dis- 
entangled in p^p —> jet7r + X. Accordingly, in the GPM 
(keeping only the leading contributions after integration 
over the initial intrinsic transverse momenta [22l | ) the nu- 
merator of the SSA for inclusive production of leading 
pions inside a large pt jet can be written as 



Ej dAa 



cPpjdzcPk^n s 
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where Ej and pj are the energy and momentum of the 
observed jet, x a ,b and fc± a ,fc = /c± a ,h(cos0 ai 6,sin0 a: 6,O) 
are the initial parton light-cone momentum fractions and 
intrinsic transverse momenta respectively, and z and fej_„. 
(k±Tv — |fc_Lu-|) are the light-cone momentum fraction and 
transverse momentum of the pion inside the jet with re- 
spect to the jet (parton c) direction of motion. 

The second line of Eq. (fTJ corresponds to the Sivers ef- 
fect, with the azimuthal modulation sin0s A , where 4>s A 
is the angle of the transverse spin vector, Sa, of hadron 
A (with mass M), relative to the jet production plane 
and fvfixa^k 2 ^) is the Sivers function. H^ b ^ cd (s,i,u) 
is the unpolarized squared hard scattering amplitude for 
the process a b — > cd, with s,t,u the usual partonic 
Mandelstam variables and D\(z, k 2 ^) is the unintegrated 
fragmentation function for parton c to fragment into a 
pion (with mass M^). The third line of Eq. (jTJ) corre- 
sponds to the Collins effect, with the azimuthal modula- 
tion sin(05 4 — 4>n)i where <j>^ is the azimuthal angle of 
the pion three-momentum around the jet thrust axis, as 
measured in the fragmenting parton helicity frame, and 
H^~ c (z, kj_ ) is the Collins function. It is convoluted with 
the unintegrated transversity distribution, hf(x a , fc^ a ), 
i.e. the distribution of transversely polarized quarks in a 
transversely polarized hadron. d^N is the partonic spin 
transfer asymmetry for the process a^b — > c^d (defined 
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its azimuthal phase (for details see 22]) 
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In close analogy with the case of SIDIS, one can define 
azimuthal moments and project out the various angular 
modulations in terms of (f>s A and (j>^ from Eq. (TT]): 
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where dAa = [da(4>s A ,4>") - da((f>s A + tt, as 
given in Eq. © and da unp = [da{<j) SA ,<j>") + da{<j>s A + 
tt,(/)^)]/2 is the unpolarized cross section. By choosing 
W((j>s A ,4>n ) = sm 4>s A one then singles out the Sivers 
contribution to , that is A*^ 1 ^ Sa , which we focus on 
in the following. Moreover, since our aim is to study the 
process dependence of the quark Sivers function, we will 
consider pion-jet production at large rapidities. Here, 
any potential sea-quark and gluon Sivers effects are ex- 
pected to be negligible, as follows from the analyses of 
SSAs in SIDIS (see 0) and in pp ->• tt + X at midra- 
pidity H [5(3, H3] and from the study carried out in [iH . 

In the GPM, the Sivers function is assumed to be uni- 
versal and taken to be the same as that probed in SIDIS; 
that is in Eq. {U, 



a.SIDIS/ ,2 
IX 1 j a* j 



(2) 



On the other hand, for the process p^p — > jet tt + X, both 
ISIs and FSIs contribute and thus in principle the Sivers 
function for the pion-jet production should be different 
from that probed in SIDIS. Following we carefully 
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FIG. 1. The Sivers asymmetry A N for the p p — > jet7r + X as a function of pjr, at fixed jet rapidity nj = 3.3, for RHIC 
energy, y/s = 500 GeV. The solid (SIDIS 1 HjJ) and dotted (SIDIS 2 curves are for the GPM calculation, and the dashed 
and dot-dashed ones for the CGI GPM calculation. The vertical dotted line corresponds to if = 0.3. 



analyze these ISIs and FSIs for all the partonic scatter- 
ing processes relevant to the azimuthal distribution of 
leading pions inside a fragmenting jet in proton-proton 
scattering. In this way the (quark) Sivers contribution 
from the CGI GPM will be (here a = q) 



In Fig. Q] we plot A™* Sj4 (pj) for 7r°' ± -jet production, as 
a function of the jet transverse momentum p^T at for- 
ward rapidity, rj 3 — 3.3, for RHIC energy, Jli_ = 500 
GeV, integrated over k± n and z (z > 0.3) [22]. The 
estimates using the two available parameterizations of 
the Sivers function in the GPM formalism are shown as 
the solid (SIDIS 1 [29]) and dotted (SIDIS 2 gU) lines, 
while the corresponding ones using CGI GPM formalism 
in Eq. ((3]) are shown as dashed and dot-dashed lines. One 
immediately sees that the results of the two approaches, 
while comparable in size, exhibit different signs. The 
opposite sign is the manifestation of ISIs and FSIs, see 
Fig. [2] for illustration. Particularly for the dominant 
channel at forward rapidity, qg — ¥ qg with the final quark 
identified with the observed jet, these ISIs/FSIs lead to 

H lT^9 ~ zIfiS fOT CGI GPM ' 

for GPM fill/ 
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in which a process- dependent Sivers function denoted as 
j±a,ab->cd ^ usec j ra ther than that from SIDIS as in the 
GPM approach 2lj, [22J. The crucial point is that the 



existence of the Sivers function in the polarized nucleon 
relies on the ISIs and FSIs between the struck parton and 
the spectators from the polarized nucleon through the 
gluon exchange. Thus by analyzing these interactions, 
one can compute the color factors Cj (Cf c ) for initial 
(final) state interactions that determine the process de- 
pendent Sivers function to be used for the corresponding 
partonic scattering ab — > cd. In the CGI GPM, the pro- 
cess dependence of the Sivers function can be shifted to 
the squared hard partonic scattering amplitude, that is 



/ 
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where all the process dependence is absorbed into the 
new hard function H a ^ cd , which is the same as in the 
single inclusive particle production This approach 

suggests a close connection with the twist-3 collinear for- 
malism [ll], [13] (see [HI for details) . 

Now we study the consequence of these ISIs and FSIs 
by comparing the predictions of the Sivers asymmetry 
for pion-jet production between GPM and CGI GPM. 



while 

qg- 



►93 





FIG. 2. Sample Feynman diagrams for the initial (a) and final 
(b) state interactions for pion-jet production, illustrated here 
for the partonic channel qg — > qg with the final q identified 
with the observed jet. 

The predictions labeled SIDIS 1 and SIDIS 2 are simi- 
lar in the intermediate pjt < 5.5 GeV region (correspond- 
ing to Feynman x, xp < 0.3), where the Sivers function 
parameterizations are constrained by present SIDIS data 
(that is at Bjorken x, xb < 0.3). This region is then 
optimal to test directly the process dependence of the 
Sivers function. This is the main goal of our analysis. 
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Moreover, the observation of a sizable A N A at large 
PjT (i-e. large xp) could be extremely useful to constrain 
the Sivers function in the large x region 22], as well as 
to test its process dependence. In this respect large Xb- 
data from SIDIS [3(3], e.g., at Jefferson Lab in the 12 GeV 
program, could be very important. 

Note that for yfs — 200 GeV the behavior of our esti- 
mates would be similar to that shown in Fig. [I] gaining 
almost a factor of 2 in size. However the range of pyr 
covered would be narrower (pjx < 6.5 GeV) and with 
xf < 0.3 now corresponding to p- 3 T < 2.2 GeV. 

As a natural extension of this work we can consider 
single inclusive jet asymmetry in proton-proton scatter- 
ing by replacing the fragmentation function D\{z,h\^) 
in Eqs. (|T]) and ©, by 5{z - l)S 2 (k ±n ). In this case 
the SSAs are described solely by the Sivers function. An 
analogous study of the Sivers contribution yields a similar 
process dependence and the results we obtain for A a ^ ^ Sa 
(not shown) look almost indistinguishable from the case 
of neutral pion-jet production (central panel of Fig. [T]). 

In summary, we have studied the azimuthal distribu- 
tion of leading pions inside a jet as well as single in- 
clusive jet production in proton-proton scattering, under 
present active investigation at RHIC. By adopting the 
TMD GPM, we have considered ISIs and FSIs leading 
to process dependence of the Sivers function. We have 
presented estimates of the Sivers asymmetry A s ^^' Sa 
for RHIC kinematics within the GPM framework with 
and without inclusion of color gauge factors. We find 
that the resulting Sivers asymmetries are comparable 
in size but appear with opposite signs. We conclude 

sin (f) s 

that the experimental observation of a sizable A N A 
in p1~p — > jet 7r + X or p^p — > jet + X can test the role 
color gauge invariance plays in the universality properties 
of the Sivers function. At the same time it could give a 
clean indication on the size of the Sivers function in the 
large x region (not covered by present SIDIS data). This 
will definitely provide new insights into our understand- 
ing of single spin asymmetries in QCD. 
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